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I. INTRODUCTION AND NOTATIONS 



In this note we adopt the space-time relativistic approach of Feynman's propagators (for bosons and fermions) 
instead of the canonical Lagrangian-Hamiltonian quantized field theory. The rationale for this choice is that for the 
development of our basic ideas the former alternative is better suited to exhibit the connection between the propagator 
of quantum mechanics and the underlying Levy processes. More precisely, the relativistic Feynman propagators are 
here linked to a dynamical theory based on a particular Levy process: a point, already discussed in a previous 
^ ' paper [l[ , which is here analyzed thoroughly with the purpose of deducing its consequences for the basic interactions 
among the fundamental constituents, namely quarks, leptons, gluons, photons etc. To this end we first recall that a 
■ Levy process is a stochastic process X(t), t > on a probability space (il, F, P) such that 

(N ; • X(0) = 0, P-q.o. 

• X(t) has independent and stationary increments: for each n and for very choice of < t\ < t2 < • • • < t n < +oo 
the increments X(tj+i) — X(tj) are independent and X(tj+\) — X(tj) = X(tj + i — X(tj); 

• X(t) is stochastically continuous: for every a > and for every s 

\unP(\X(t)-X( S )\>a)=0. 

b : 

5~ , To simplify the notation we will consider in the following only one-dimensional Levy processes (an n-dimcnsional 
extension, however, would not be a very difficult task): it is well known [2-4] that all its laws are infinitely divisible, 
but we will be mainly interested in the non stable (and in particular non Gaussian) case. In other words the 
characteristic function of the process At-increment is [ip(u)] At ^ r where <p is infinitely divisible, but not stable 1 , and 
r is a time scale parameter. The transition probability density p(2|l) = p{x2, t2\xi, t\) of a particle moving from the 
space-time point 1 to 2 then is 

1 r + °° 

p(2\l) = — du[(p(u)] ( -*'- tl V T e- iu( - x '- xl ) (1) 
2tt J-oo 



'Electronic address: cufaro@ba.infn.it 

^Electronic address: pustcrla@pd.infn.it 

A law ip is said to be infinitely divisible if for every n it exists a characteristic function ip n such that ip = </p"; on the other hand it is 
said to be stable when for every c > it is always possible to find a > and ieR such that e' bu ip(au) = [ip(u)] c . Every stable law is 
also infinitely divisible. 
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In analogy with the non rclativistic Wiener case, we then obtain for the motion of a free particle the Feynman 
propagator /C(2| 1) = IC(x2,t2\xi,ti) as 

£(2|1) = — / du [^(u)] i ( i2 -* 1 )/ T e- i " (x2 - a:i) (2) 
and the corresponding wave function evolution is 

/+oo 
dx' JC(x,t\x',t')ip(x',t'). (3) 
-oo 

From ([2]) and ([3]) we easily obtain [l[ 



id t ip = — vidx)^ 

T 

where r) = log ip and rj(d x ) is a pseudodifferential operator with symbol r)(u) d, [B-Q which plays the role of the 
generator of the semigroup Tt = e tr, ^ d ^/ r operating on a Banach space of measurable, bounded functions 0, 0-0] ■ 

It is very well known 0, on the other hand, that ip represents an infinitely divisible law if and only if rj(u) = 
log <p(u) satisfies the Levy— Khintchin formula 

B 2 u 2 f 

r)(u) = iju - h / [e mx - 1 - iux I[_ 1A] (x)] v(dx) (4) 



where j,/3 6 M, I[_i ^(x) is the indicator of [—1, 1], and v(dx) is a Levy measure, namely a measure on K such that 
u({0}) = and 

(a; 2 A 1) v(dx) < +oo. 

In the most common cases of centered and symmetric laws the equation (j4]) simplifies in 

B 2 u 2 f 

rj(u) = h / (cosux — 1) v(dx) (5) 



and r){u) becomes even and real. As a consequence the corresponding operator r](d x ) is self-adjoint and acts on 
propagators and wave functions according to the Levy-Schrodinger integro-diffcrential equation 

id t ^(x,t) = --r](d x )^{x,t) = -^-dl^(x,t)-- [ [i/>(x + y,t)-1>(x,t)]v(dy) (6) 

The integral term accounts for the jumps in the trajectories of the underlying stochastic process, while an action a with 
B 2 = ar/m weights the usual differential term of the Schrodinger equation. For (B — a pure jump Levy-Schrodinger 
equation is obtained 

id t i>{x, t) = -- [ [i>{x + y,t)- t)\ v{dy). (7) 

II. STATIONARY SOLUTIONS FOR THE FREE PARTICLE 

The free equation ([B]) admits simple stationary solutions: if we consider 

Mx, t) = e~ lEot / a cj)(x), a = 

T 

we then have 

cv 2 cv r 

Eo4>(x) = -=-f{x) - - / [4>{x + y) - <t>{x)} v(dy), (8) 
2m r J M 
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and for a plane wave <fi(x) = e ux from ([5]) - namely with a symmetric Levy noise - we get 



2„,2 



+ / (e ±<u » - 1) Kdy) 



2, ,2 



/3 2 U 



+ / (cos uy — I) is (dy) 



(x) = --rj(u)<j){x) 



which is satisfied when Eq = —ctr](u)/T. Hence, by taking p = au as a momentum variable, we obtain [l[ the relevant 
equation 



E, 



o 



t \a/ 



which connects the kinetic energy of a forceless particle to the logarithmic characteristic of a Levy process. 



(9) 



III. RELATIVISTIC QUANTUM MECHANICS 

Let us take now in particular the non stable law 



rj(u) = 1 - \j\ + a 2 u 2 . (10) 
With the following identification of the parameters 

k H 2 h n 

a = n, — = mc , a — — , p — nu. 

t mc 

we are led to the formula 

Eo = —mc 2 r\ ( ^ ) = E — mc 2 = \J m 2 c i + p 2 c 2 — mc 2 (11) 



which is the well-known relativistic kinetic energy for a particle of mass m. The Schrodingcr equation of a relativistic 
free-particle is then easily obtained from (fTl~T) by reinterpreting as usual E and p respectively as the operators ihd t 
and —ihd x 



ihd t ip{x,t) = v / m 2 c 4 - h 2 c 2 d 2 x ip(x,t), (12) 

It is easy to check that this derives also from ((6]) after absorbing the mass energy term —mc 2 of (fTTj) into a phase 
factor e lmc t ' fi . Remark that in three dimensions (fT2"j) would read 



iHd t ip(x,t) = yjm 2 c± - h 2 c 2 V 2 %k{x,t) (13) 

It has been shown 0, Q that the Levy process behind the equations (fT2j) and (fT3|) is a pure jump process 0, H[ with 
an absolutely continuous Levy measure v(dx) = W(x)dx such that 

W( X )='K 1 (m='K 1 (^\ X \) (14) 
7r|a;| \ a J tt\x\ \ n > 

(K v are the modified Bessel functions 0), that in three dimensions becomes 

, \ 1 /W\ mc I mc , \ 

- K 2 (U. j = K 2 (- W ) (15) 

while from ||7J) the equation (fT2"|) takes the form 

u + \ 2 f i/j(x + y,t)-ip(x,t) (mc A 

ihd t ^(x, t) = -mc / n Ai I -r- \y\) dy (16) 

and in three dimensions is 

■%x it +\ 2 /* ip(x + y,t) -ip(x,t) mc (mc \ 3 

iHd t ip(x,t) = -mc A \ — X 2 — y cfy 17) 
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From the equation (|T3|) ~ by means of well known standard procedures [10| - one also derives (always for the free 
particle) the Klein-Gordon and Dirac equations in three dimensions both for scalar, and for spinor wave functions ip, 
namely respectively 

(□-^)tf = 0, (18) 

(i7^ - ™) V = 0. (19) 

The corresponding Klein-Gordon and Dirac propagators in their turn satisfy the inhomogeneous equations (here with 
h = c = l) 

(□a - m 2 ) K KG (2\l) = <5 4 (2|1) (20) 
(ry^ -m)K D (2|1) = z<5 4 (2|l) (21) 

with 5 4 (2|1) = 5(t2 — ti)8 3 (x2 — Xi). Let us finally remark that these relativistic quantum wave equations have been 
recently of particular interest also in the field of quantum optical phenomena and of quantum information. 



IV. INFINITE DIVISIBILITY-PRESERVING MODIFICATIONS 



From the relativistic kinetic energy Eq of a point particle of rest mass m 



E0{P) --/l + (22) 



with the identifications 



we obtain 



E P 
u = 



mc 2 ' amc 



V {u) = 1 - yT 



a 2 u 2 



which of course coincides with the equation (|10[) extended to three dimensions. We take now a class of transformations 
of 77(11) characterized by the fact that they preserve the infinite divisibility, while producing changes in the forceless 
particle equations of motion with respect to the usual ones (|18[) and (|19j). To this end we modify the energy-momentum 
formula as follows 




E(p) = mc\ ! + -£— + / (23) 



where / is a - possibly small - dimensionless, smooth function of the relativistic scalar p 2 /m 2 c 2 . Of course this 
modification entails that p 2 no longer coincides with m 2 c 2 since the standard energy- momentum relation is now 
changed into 

2 E 2 2 22,22c(v 2 \ r n ,\ 

V =— P =m c + m c / —j-g . (24) 



As we will see in the following, this also implies that the mass no longer is m: it will take instead one or more values 
depending on the choice of /. As a matter of fact, it could appear to be preposterous to introduce a function / of 
an argument which after all is a constant (albeit different from 1). However we will show that this artifice will lend 
us the possibility of having both a mass spectrum, and a new wave equation when - in the next section - we will 
quantize our classical relations. Moreover it will be argued in the following that we will choose / in such a way that 
the corresponding modified logarithmic characteristic 77 will remain infinitely divisible: a feature that is instrumental 
to keep a viable connection to a suitable underlying Levy process. 

To see that, we first remark that ([24]) defines the total particle energy E in an implicit form. To find it explicitly 
we first rewrite (|24|) in a dimensionless form as 

J?_ = 1 + f f P 2 
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and then, by taking g(x) — x — f(x), we just observe that the former equation requires that x be solution of g{x) = 1, 
namely 



fJ 



P 



P 



f 



P 



1. 



Note that / and g should be considered universal functions, and that the following conditions hold 

/(1) = = 1 

If then represents one of the (possibly many) solutions of this equation, we could write 



P- 



g- 1 (l) = x(g) 



9=1 



so that we have 



p 2 = ^-p 2 = m 2 c 2 g-\l) 



which can be interpreted as a simple mass re-scaling, from m to one of the (possibly many) values M = m^J g -1 (l). 
The new hamiltonian then is 



E(p) = v / m 2 c' i g- 1 (l) +p 2 c 2 = Mc 2 \ 1 



P 



M 2 c 2 



(25) 



and its kinetic part (by applying the same re-scaling also to the subtracted rest mass term) is 

Eo(p) = E(p) - mc 2 ^g~Hl) = Mc^l + j^ - Mc 2 . 

Hence the main consequence of our modification consists of a re-scaling of the mass value (to — > M) at a purely 
classical level. This fact is apparently welcome, because its straightforward consequence is that the new associated 
logarithmic characteristic 77 is trivially again infinitely divisible, and hence still produces acceptable Levy processes. 
But there is more: since g~ l {l) can take several different real and positive values, by means of our modification 
we have introduced a mass spectrum: in the rest frame of the particle we indeed have now 



M = E cm /c 2 = m^g- l (l) 



(26) 



V. QUANTUM EQUATIONS OF MOTION 



From the modified energy formula one derives a relativistic Schrodinger equation (for instance by means of the 
formal substitutions E — > ihd t and p — s- —ihV): 



iHd t ip(x,t) = mc 2 J I - -^V 2 + f (-^2) Hx,t) (27) 

where the square root pseudo-differential operator satisfies the constraints exposed in the Section [U From (|27[) one 
easily obtains in the usual manner (from here on K = c = 1) 



□ - m / I □ I - m 

TO 



□2 - m 2 f 



\ TO' 



□ 2 " 



^ = 0, 

K KG {2\1) = 5\2\l) 



and, by standard methods [10j . the modified Dirac spinor equations 



□ 



ip = 



i7 ^-T +/ (i D2 ) 



K D (2\1) = z<5 4 (2|l) 



(28) 
(29) 

(30) 
(31) 
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In the momentum space (with Fourier transforms in four dimensions) these equations become much simpler: more 
precisely we have 

1Ckg(p 2 ) = -o 2 h . ff 2 / 2 M i ( 32 ) 

p z - m z [1 + j \p z i 'm z )\ + it 

JC D {p 2 ) = 1 (33) 

7 M £V» — m v 1 + f(p 2 / m2 ) + i e 

We notice that /C_d(2|1) is in our case simply related to the /C_r:g(2|1) (like in the usual case) as 

£d(2|1) = i [i?2 + my/1 + f(D 2 /m 2 ) ] 1C KG (2\1) (34) 

VI. PHENOMENOLOGY: QUARK AND LEPTON MASSES 

The equations (f2"§)) and (f3"Tj) generalize the well known propagator equations ([2H)l and (|21[) deriving from QED and 
QCD at zero order (in absence of interaction terms). The Standard Model (SM) 2 SU C (3) x 5J7l(2) x U(l) treats both 
strong, and electro-weak interactions: within this scheme the modified 77(11) leads to new interesting consequences. 
We begin by considering the Feynman rules in perturbation theory in presence of the modified zero order propagator 
for both spin i (quarks and leptons) and spin 1 (gluons, vector weak interacting Bosons). The amplitude A for a 
fermion that propagates from vertex X to vertex Y, if expanded, looks as follows: A = A^ + A^ 1 ' + A^ + . . . The 
lowest order is 

A<® = Y 1 X. (35) 

7 M £V - myl + f(p 2 /m 2 ) + ie 

It is then possible that the Fermion emits and reabsorbs a virtual vector boson 3 from X to Y: 

7 m 1 



A« = AnglY J d 4 k ■ 

X (36) 



YPp ~ my/l + f{p 2 /m 2 ) (jp - k) 2 

1 7m 



k"^ v - my/1 + f(k 2 /m 2 ) + ie jPp p - m^/l + f(p 2 /m 2 ) 

We choose now f{x) in such a way that it makes finite the integral 

f d A k 1 

G ~ ^ J Yk p - my/T+ f (k 2 /m 2 ) +i~e {p - k) 2 7 " (3? ^ 

One may notice that f(x) behaves as a smooth cut-off in a procedure of regularization at each order in QCD (and 
QED). The integral C is an invariant of the form C = A(p 2 )%f + B(p 2 ) and its integrand is also present as a factor in 
higher order terms, thus producing convergence. 



2 For future developments we recall that the Lagrangian density of QCD is, up to gauge fixing terms: 

1 

where = dpA® — d„A° + gefabcA^A^, and the insertion of interaction terms is done with the minimal interaction by substituting 
the simple derivative dp with the covariant one Dp where we have respectively for QED and QCD 

Dp = dp- ieAp 
(Dp) l3 = S^dp-igsT^Al. 

Here g s is the QCD coupling constant, and f a bc are the SU(3) color matrices and structure constants respectively, and A^ the eight 
Yang-Mills gluon fields; ip? are the Dirac 4-spinors associated with each quark field of color i and flavor q. 

3 The presence of gammas in the numerator of formulae II36H and J 371) is typical of QED. More elaborated numerators may be present in 
non-abelian theories (in particular QCD); however they appear totally unessential for our subsequent developments and purposes. 



Cufaro Petroni and PusterlA: Masses and Levy-Schrddinger equation 



The search of poles of the fermion propagators can be done in the following way QJl : one considers the contributions 
of the perturbative expansion of the amplitude A(p 2 ) (here we always understand / = f(p 2 /m 2 )): 



A(p 2 ) = Y 



1 



1 



-,C- 



1 



my/1 + f y — m VI + / y — my/1 + f 
1 „ 1 



--C 



c- 



-ms/T+J y-7WTT7 tf-mVTTJ 



and using the formula 



1 111111 
-B— + -B—B- 



A-B A A A A A A 
one obtains the approximate expression 

1 „ „ 1 



A~Y- 



-X = Y- 



-^X 



X 



my/T+f -C / - my/TTJ — A%f - B 

and looks for possible poles which - after rationalizing equation (|40[) - are solutions of the equation 

= 



1-A(p 2 ) p 2 - my/l + f(p 2 /m 2 ) + B(p 2 ) 



(38) 



(39) 



(40) 



(41) 



A. Hypothesis for an approximate evaluation of the mass spectrum 

Let us now focus our attention on QCD. We may consider the approximation A(p 2 ) ~ A(m 2 ) and B(p 2 ) ~ B(m 2 ) 
which follows from the assumption 



A < 1 and B < my/i + f(p 2 /m 2 ) 



We then obtain the equation 



my/1 + f{p 2 /m 2 ) + B(m 2 ) 
1 -A{m 2 ) 



(42) 



(43) 



where the experimental masses m exp are still represented in an implicit form. However note that in the limit g s — > 0, 
A and B — > one achieves the equation 



p 2 = m 2 [l + f(p 2 /m 2 )] 



(44) 



which coincides with the classical equation (|24j) . At this point we notice that the simplest choice of f(x) that makes 
the integral C finite (integrand convergent) is a polynomial of third degree: 



f(x) = A + Aix + X 2 x 2 + X 3 x 3 

The equation (|44| then becomes 

x — 1 — f(x) = — Xs(x — l)(x — x+)(x — X-) = 



with 



where 



/(I) = o 



x± = — (-A 2 - A 3 ± VA 



Aq — — Ai — A2 — A3 



A = (A 2 - A 3 ) 2 - 4AxA 3 - 4X 2 3 + 4A 3 



(45) 



(46) 



(47) 



(48) 
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and we finally get 



A2 
A 3 
Ai-1 
A3 



x + + X- + x + x 



(1 +X+ +X-) 



(49) 



(50) 



From the previous formulae one achieves the following interesting result: the convergence of C determines a possible 
phcnomcnological function f(x) that produces a mass spectrum of three fcrmion particles (quarks in QCD). 

The connections with the possible experimental physical masses are Mi = m, M2 = m^Jx+, M3 = m^Jx- . If the 
three poles in the free (zero order) propagator are real and positive (with proper residues) , with appropriate values of 
the A's, they allow the interpretation of physical basic masses of fermions (quark or leptons) belonging to the three 
different families of the Standard Model. To be more specific we get two different propagators for quarks, one with 
charge — | (d,s,b quarks) and another with charge +| (u, c, t quarks). Similarly for charged leptons (charge —1 and 
spin i) we get one propagator. 



We have proposed a modification of the classical relativistic hamiltonian that allows the presence of several masses 
without changing its basic structure. This modification does not affect the infinite divisibility of the laws that 
are at the basis of the correspondence between stochastic processes and Levy-Schrodinger equations. However we 
discovered that the mentioned modification suggests a reformulation of the relativistic equations for wave functions 
and propagators in such a way that a suitable choice of the background noise produces a convergence in the perturbative 
contributions. To this purpose we remarked that a modification - with respect to the one given by equation 
- of the logarithmic characteristic f](u), by the insertion of a cut-off f(x), allows to proceed to rcgularization first, 
and then to renormalization of the two-point function of QCD. There are three parameters in our phenomcnological 
f(x) which is a third degree polynomial; the latter appears as the simplest choice that produces convergence in the 
integrals representing high order contributions to the fermion and boson propagators. Such parameters create three 
different poles in the zero-order propagators and allow the interpretation of a physical system with three different 
masses under precise constraints on f(x). The masses might be related to the three families of the Standard Model. 
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